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1 Introduction 



A fundamental problem in game theory is that of finding a Nash equilibrium of a 
bimatrix game, that is, a two-player game in strategic form. This is achieved by the 
classical pivoting algorithm by Lemke and Howson (1964). Shapley (1974) intro- 
duced the concept of an index of a Nash equilibrium, and showed that the endpoints 
of every path computed by the Lemke-Howson algorithm have opposite index. As 
a consequence, any nondegenerate game has an equal number of equilibria of posi- 
tive and negative index, if one includes an "artificial equilibrium" (of, by convention, 
negative index) that is not a Nash equilibrium. The Lemke-Howson algorithm is one 
motivating example for the complexity class PPAD defined by Papadimitriou (1994). 

Euler complexes, introduced by Edmonds (2009), provide a more recent abstract 
framework for the Lemke-Howson algorithm. A J-dimensional Euler complex (or 
"J-oik") over a finite set of nodes is a multiset of ^^-element sets called rooms so that 
any set of <i — 1 nodes is contained in an even number of rooms; if these are always 
zero or two rooms, this is the familiar concept of an (abstract simplicial) manifold. 
For a family of oiks over the same node set V , Edmonds (2009) showed that there 
is an even number of room partitions of V . A special case is a family of two oiks 
of possibly different dimension corresponding to the two players in a bimatrix game. 
Then room partitions are equilibria, and the Lemke-Howson algorithm is a special 
case of the "exchange algorithm" used to show that there is an even number of room 
partitions. In another special case, all oiks in the family are the same 2-oik, which is 
an Euler graph with edges as rooms and perfect matchings as room partitions. 

This paper presents three main contributions in this context. First, we define a 
unifying formalism called pivoting systems that describes "complementary pivoting 
with direction" in a canonical manner. Second, using this formalism, we extend the 
concept of orientation to oiks and show that room partitions at the two ends of a 
pivoting path have opposite sign, provided the underlying oik is oriented. Third, 
a room partition in the special case of an oriented 2-oik corresponds to a perfect 
matching of a graph with an Eulerian orientation. Here we give a polynomial-time 
algorithm to find another perfect matching that has opposite sign (the complementary 
pivoting algorithm that achieves this may take exponential time). 

Our concept of a pivoting system (see Definition 2) has the following features. 
A pivoting operation switches back and forth between two states by changing one 
component of an m-tuple representing each state. The system is oriented if these 
two states have always opposite orientation. An example of a state is a vertex of a 
simple polytope, represented by the m facets it lies on, and oriented by the sign of 
the determinant of the normal vectors of these facets. Labels define "complementary 
pivoting" to connect completely labeled states by paths of intermediate almost com- 
pletely labeled states. In an oriented system, the endpoints of a path have opposite 
sign, and the direction on the path can be locally identified by two signs associated 
with an almost completely labeled state (see Theorem 3). We think this approach 
captures "directed complementary pivoting" in the most natural way. 
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Our framework is similar to that of Lemke and Grotzinger (1976), except that we 
distinguish between a state and its representation to capture more general cases such 
as room partitions. Todd (1974; 1976) considered abstract complementary pivoting 
with "primoids" and "duoids", where the latter are the same as the "oiks" by Edmonds 
(2009). However, the room partitions implicit in Todd's approach are limited to two 
rooms. Eaves and Scarf (1976) generalize the algorithms by Lemke and Howson 
(1964), Lemke (1965), Scarf (1967), and others, as following paths defined by a 
piecewise linear map. They extend the index theory of Shapley (1974) to these paths; 
our use of determinants (Proposition 4) to obtain the orientation is much simpler. 

Orientation in oiks, which we study in Section 3, seems to be a new concept (see 
Definition 6), which extends the known definition for manifolds (e.g., Hilton and 
Wylie, 1967). The "sum" of oriented oiks is then again oriented. This implies that 
ordered room partitions come in pairs of opposite sign (Theorem 10). Keeping the 
order of rooms is necessary for oiks of odd dimension, for example when the rooms 
are (abstract) triangles, as discussed for the example of an octahedron in Figure 1. 
For even dimension, the order of the rooms in a room partition can be disregarded (as 
indeed for the complementary pivoting path), while keeping the property of oppo- 
sitely signed room partitions at the end of the path (Theorem 1 1 , proved via pivoting 
systems). Instead of labels, one can use room partitions with a "Sperner oik" (Ed- 
monds, Gaubert, and Gurvich, 2010); the connection is very close, as we discuss in 
Appendix A. 

Section 4 is concerned with signed perfect matchings, which for a graph with an 
Eulerian orientation are exactly the signed unordered room partitions in an oriented 
2 -oik. The sign of a matching is the parity of the permutation of the nodes of the graph 
when writing down the matched edges as they are oriented. Signs of matchings have 
been extensively studied in the context of Pfaffian orientations, which are orientations 
of a graph so that all matchings have the same sign. Eulerian orientations are not 
Pfaffian because they have equal numbers of matchings of either sign. This follows 
also from the fact that the skew symmetric incidence matrix of the graph is singular 
in that case, so that its determinant and hence the Pfaffian of the matrix is zero, as 
mentioned in our exposition of these results at the beginning of Section 4. 

The question whether a general orientation of a graph is Pfaffian is polynomial- 
time equivalent to deciding if the graph has a Pfaffian orientation (see Vazirani and 
Yannakakis, 1989, and Thomas, 2006). For bipartite graphs, this problem is equiva- 
lent to finding an even-length cycle in a digraph, which was long open and shown to 
be polynomial by Robertson, Seymour, and Thomas (1999). For general graphs, its 
complexity is a notoriously difficult open question. 

For a graph with an Eulerian orientation and a given matching, there is another 
matching of opposite sign. Its existence is guaranteed by the complementary pivoting 
algorithm, which, however, may take exponential time. In Theorem 12 we give an 
algorithm to find such an oppositely signed matching in polynomial time. It makes 
essential use of the Euler property, because other approaches seem to lead to the 
difficulties associated with Pfaffian orientations in general graphs. Merschen (2012, 
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Theorem 5.3) has shown how to find in polynomial time an oppositely signed match- 
ing in a planar Euler graph, and his method can be adapted to graphs that, like planar 
graphs, are known to have a Pfaffian orientation. Our algorithm for general Euler 
graphs is surprisingly simple and can be implemented in near-linear time in the num- 
ber of edges of the graph; this implementation is presented in detail in Appendix B. 

For bipartite Euler graphs, it is even simpler to find a matching of opposite sign, 
using the complementary pivoting algorithm. We extend this linear-time algorithm to 
bipartite graphs that are oriented so that no node is a source or sink (Proposition 13). 

In general, complementary pivoting for perfect matchings in graphs with an Eu- 
lerian orientation may take exponential time. As discussed at the end of Section 4, 
Casetti, Merschen, and von Stengel (2010) and Merschen (2012) have shown that one 
can use the exponentially long "Lemke paths" of Morris (1994) for this purpose. Ed- 
monds and Sanita (2010) describe 3-manifolds where finding a second room partition 
may take exponential time. 

Section 5 concludes with a discussion of the computational complexity of pivot- 
ing systems. 

2 Labeled oriented pivoting systems 

This section presents "pivoting systems" as a unifying formalism for path-following 
with "complementary pivoting". This abstracts from the "Lemke paths" on labeled 
polytopes as applied to bimatrix games and linear complementarity problems, which 
we describe first. In an oriented pivoting system, the endpoints of complementary 
pivoting paths have opposite sign, for which Theorem 3 provides a canonical proof. 
Polytopes are oriented via signs of determinants (Proposition 4). 

We use the following notation. Let [fc] = { 1 , . . . , /c} for any positive integer k. The 
transpose of a matrix B is . All vectors are column vectors. The zero vector is 0, 
the vector of all ones is 1, their dimension depending on the context. Inequalities 
like X > between two vectors hold for all components. A unit vector e^ has its ^th 
component equal to one and all other components equal to zero. A permutation jt 
of [m] has parity ( — 1)^ if k is the number of its inversions, that is, pairs /, j so that 
/ < j but n{i) > n{j), and the permutation is also called even or odd when k is even 
or odd, respectively. 

A polyhedron P is the intersection of n half spaces in R*", 

P = {xeM."' \a]x<bj, j e[n]} (1) 

with vectors aj in R"* and reals bj. A labeling function /:[«]—> [m] assigns a label 
to each inequality in (1), and .x: in P is said to have label l{j) when the jth inequality 
is binding, that is, ajx = bj, for any j in [n] . 

We normally look at the "nondegenerate" case where binding inequalities define 
facets of a polytope, and no more than m inequalities are ever binding. That is, 
we assume P is a polytope (that is, bounded) which is simple (every vertex lies on 
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exactly m facets) and that none of the inequaUties can be omitted without changing 
the polytope, so for every j in [n] the 7th binding inequaUty defines a facet Fj given 
by 

Fj = {xeP\ ajx = bj} (2) 

(for notions on polytopes see Ziegler, 1995). Then facet Fj has label for j in [n] , 
and we call P a labeled polytope. A vertex of P is completely labeled or CL if the m 
facets it lies on have together all labels in [m] . 

A linear complementarity problem (LCP) with an m x m matrix M and m -vector 
q is the problem of finding a CL point z of the polyhedrom 

P = {z e M"* I -z < 0, -Mz < ^} (3) 

whose 2m inequalities have labels l,...,m,l,...,m. A CL point z is also called 
"complementary" because for each / in [m] either Zi = or {q — Mz)i = 0. Lemke 
(1965) described a path-following method of "complementary pivoting" to solve cer- 
tain LCPs. If P in (3) is a simple polytope, then CL vertices are the unique endpoints 
of these "Lemke paths" (see also Morris, 1994) and hence there is an even number of 
them; we prove this in more general form in Theorem 3 below. 

Lemke paths also correspond to the paths computed by the algorithm by Lemke 
and Howson (1964) that finds one Nash equilibrium of a bimatrix game (see von 
Stengel, 2002, for a survey). Suppose the polytope P has the form 

P={xeR'"\-x<0,Cx<l} (4) 

for some {n — m)xm matrix C, and that each of the first m inequalities > has 
label i in [m]. Then is a completely labeled vertex. If P in (1) has a completely 
labeled vertex, it is easy to see that it can be brought into the form (4) by a suitable 
affine transformation that maps that vertex to 0. If C is a square matrix, then the CL 
vertices x of P other than correspond to symmetric Nash equilibria (jc',jc') of the 
symmetric game with payoff matrices (C, C^), where x' = x/l^x; in turn, symmetric 
equilibria of symmetric games encode Nash equilibria of arbitrary bimatrix games 
(see, e.g., Savani and von Stengel, 2006, also for a description of the Lemke-Howson 
method in this context). 

For a polytope with a general matrix C in (4), the following proposition shows 
that its CL vertices correspond to Nash equilibria of the "unit- vector game" (A, C^) . 
The unit vectors that form the columns of A encode the labels of P. 

Proposition 1 Suppose that in (4) each inequality xi > has label i, and the last 
n — m inequalities Cx<l have labels l{m + j) for 7 e [n — m]. Then x is a CL vertex 
ofP — {0} if and only if for some y the pair [x/ {\^x)^y) is a Nash equilibrium of the 
mx{n-m) game (A,C'^) where A = • • •e/(n)]. 

Proposition 1 is not hard to prove (see also Balthasar, 2009, Lemma 4.10, for a 
dual version). The special case when A is the identity matrix describes an "imita- 
tion game" whose equilibria correspond to the symmetric equilibria of the symmetric 
game (C,C^) (McLennan and Tourky, 2010). 



5 



In the remainder of this section we describe a general framework which can be 
applied to the "Lemke paths" on labeled polytopes and has other applications that we 
give later. We first give a narrative description and then a more formal development. 

Lemke paths traverse a sequence of adjacent vertices of a poly tope by pivoting 
steps just like the simplex method for linear programming, but with a different "com- 
plementary" pivoting rule. We use the general term states, which form a set S, instead 
of polytope vertices. Each state s is represented by an m -tuple 

r{s)^{si,...,Sm) (5) 

of nodes Si from a given set V . For a polytope as in (1), the set of nodes V is the set 
[n] that numbers its facets, and a state is a vertex of P represented by the m facets it 
lies on. 

The pivoting operation / takes a state s and replaces its /th component si of its 
representation in (5) by another element m of V to get a new m-tuple which we denote 

by {r{s) I i^u), 

((*!,.. I /-).m) = (51,. ..,5/_l,M,5,+l,. ..,*„). (6) 

We denote the resulting new state with this representation hy t = f{s, i) . The pivoting 
step is simply reversed by 5 = i) . (We will later refine this by allowing r{t) to be 
a permutation of r{s).) In the polytope, s and t are adjacent vertices that agree in all 
binding inequalities except for the /th one. 

Each node w in V has a label l{u) given by a labeling function / : V — > [m] . The 
path-following argument has as endpoints of the paths completely labeled (CL) states 
s where, given (5), {l{si) | ? G [m]} = [m]. In addition, it considers states s that are 
ACL or almost completely labeled defined by the condition {l{si) | i G [m]} = [m] — 
{w}, where w is called the missing label and the unique k so that k = l{si) — l{sj) 
for / 7^ j is called the duplicate label. 

"Complementary pivoting" means the following: Start from a CL state s and 
allow a specific label w to be missing, where l{si) = w. Pivot to the state t = /(5, /). 
Then if the new node u in (6) has label l{u) = w, then t is CL and the path ends. 
Otherwise, l{u) is duplicate, with l{sj) — l{u) for j ^ i, so that the next state is 
obtained by pivoting to f{t,j), and the process is repeated. This defines a unique 
path that starts with a CL state, follows a sequence of ACL states, all of which have 
missing label w, and ends with another CL state. The path cannot meet itself because 
the pivoting function is invertible, so the process terminates. 

We also want to give a direction to the pivoting path. For this purpose, a CL state 
will get a sign, either +1 or — 1 , so that the two CL states at the ends of the path have 
opposite sign. This sign is the product of two such numbers (again either +1 or — 1), 
namely the orientation o{s) of the state s when represented as r{s) = (^i, . . . ,5^), 
and the parity of the permutation % of [m] when writing down the nodes s\,...,Sm 
in ascending order of their labels. In the polytope setting, the orientation of a vertex 
is the sign of the determinant of the normal vectors aj of the facets Fj that contain 
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that vertex, see (12) below. The important abstract property is that pivoting from 
{s\,...,Sm) to ((51, . . .,5m) I Si — )■ u) changes the orientation, stated for polytopes in 
Proposition 4 below. 

In order to motivate the following definition, we first give a very simple example 
of a pivoting path with only one ACL state apart from its two CL states at its ends. 
Consider V = {ai,a2,a3,&i,&2} with labels /(ai) = l{bi) = 1, l{a2) = /(&2) = 2, 
/(as) = 3, and three states s^^s^^s^ with r{s^) = {ai,a2,ai,), r{s^) — (Z72,«2,«3)> 
r{s^) = {b2,b\,a^). Assume that f{s^, 1) = and f{s^ ,2) = . Then starting from 
the CL state and missing label 1 pivots to .v^ (by replacing a\ with hi), which is 
an ACL state with duplicate label 2 in the two positions 1 and 2. The next com- 
plementary pivoting step pivots from to (by replacing a2 with b\), where is 
CL and the path ends. The three states have the following orientations: o{s^) = 1, 
(y{s^) = — 1, (y{s^) = 1, which alternate as one state is obtained from the next by 
pivoting. Here, the two CL states 5^ and have the same orientation. They obtain 
their sign by writing their nodes in ascending order of their labels: This is already the 
case for r(s°), but in r{s^) the permutation 2, 1,3 of the labels is odd, so the sign of 
5^ becomes —1, which is indeed opposite to the sign of s^. 

In this example, we have chosen the representations of the states 5^, 5^5^ in such 
a way that the required pivoting steps can indeed be performed by exchanging a node 
at a fixed position; however, this may not be clear in advance: another representation 
of the three states might be (01,02,^3), (a2,«3,^2), {ci?„bi,b2). In this case, we 
still allow pivoting from 5° to by going from (01,02,03) to (&2,«2,«3) but with a 
subsequent, known permutation K to obtain the representation {a2,a-i^b2) of s^; for 
the orientation of the states, we have to take the parity of it into account. 

Definition 2 A pivoting system is given by (5, V, m, r, /) with a finite set S of states, a 
finite set V of nodes, a positive integer m, a representation function r : 5 — )■ V", and a 
pivoting function f : S x[m] ^ S. For a permutation 7t of [m] and r{t) = {ti,...,tm), 
let 

Then for each t — f{s, i) , there is a permutation 7t of [m] so that r^(t) — {r{s) | / ^ a) 
for some u in V, and f{t, 7c{i)) = s. The pivoting system is oriented if each state s 
has an orientation o{s), where (T:^— >{— 1,1}, so that 

C7(r) = -o{s) ■ parity (Tt) (8) 

whenever t = f{s, i) with 7t as above. 

The following simple example illustrates the use of the permutation 7t in Def- 
inition 2, which depends on /(s, /) and is part of the pivoting system. Suppose 

r{s) = (51,52,^3) = (1,2,3) and r{t) = (?i,?2,?3) = (2,3,4), where /(.y, 1) by re- 
placing 51 with 4. Then r'^it) = (?;r(l),^;r(2),^;r(3)) = ((*1,'^2,'V3) I 1 ^4) = (4,2,3), 
so ;r(l) = 3, k{2) — 1, ;r(3) = 2, that is, n says that sj becomes t^^f^j-^ except for the 
"pivot element" Si. Pivoting "back" gives s = f{t, :;r(l)) = f{t, 3). 
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It is important to note that the pivot operation / operates on states s, giving a 
new state t = f{sj), where / refers to the ith component s,- of the representation 
r{s) = {s\,. . . ,s,n). However, there may be different states s and s' with the same 
representation r{s) = r{s'), as we will see in later examples; otherwise, we could just 
take 5 as a subset of V" and dispense with r. This is one distinction to the formal 
approaches of Lemke and Grotzinger (1976) and Todd (1976), who also assume that 
the nodes ^i,...,*^ are distinct, which we do not require either. Furthermore, we 
do not give signs to the two equivalence classes of even and odd permutations of 
(si,. . .,sj), as Hilton and Wylie (1967) or Todd (1976), but instead consider unique 
representations r{s) , and build a single permutation n into each pivoting step. 

The pivoting system (S, V,m,r, f) is labeled if there is a labeling function / : V — > 
[m]. For (51, ...,5m) where Si e V for / in [m], let /(^i,...,*^) = (/(si), ... ,/(*„)), 
and consider this m-tuple as a permutation of [m] if l{si) ^ l{sj) whenever / j. If 
the pivoting system is oriented, then the sign of a CL state s is defined as 

sign(s) = (7(s)-parity(/(r(s))). (9) 

For an ACL state s, we define two opposite signs as follows: consider the positions 

/, J of the duplicate label in r{s) = {si,. . . ,s,n), that is, l{si) = l{sj) with i ^ j, and 
missing label w. Replacing l{si) with w in l{r{s)) then defines a permutation of [m], 
denoted by {l{r{s)) | ?" — > w), which has opposite parity to {l{r{s)) \ j —?■ w) because 
that permutation is obtained by switching the labels w and l{sj) in positions i and j. 
Let 

sign(5, /) = o{s) • parity(/(r(5)) | / — >■ w), (10) 

so 

sign(*,7) = (7(s) •parity(/(r(s)) \j^w) = -sign(s,j). (11) 

This is the basic observation, together with the sign- switching of a pivoting step stated 
in (8), to show that complementary pivoting paths in an oriented pivoting system have 
a direction. This direction (say from negatively to positively signed CL end-state) is 
also locally recognized for any ACL state on the path, as stated in the following 
theorem. Hence, for a fixed missing label w, the endpoints of the paths define pairs 
of CL states of opposite sign. The pairing may depend on w, but the sign of each CL 
state does not. 

Theorem 3 Let {S,V,m,r,f) be a pivoting system with a labeling function / : V — > 

[m], and fix we [m]. 

(a) The CL states and ACL states with missing label w are connected by comple- 
mentary pivoting steps and form a set of paths and cycles, with the CL states as 
endpoints of the paths. The number ofCL states is even. 

(b) Suppose the system is oriented. Then the two CL states at the end of a path have 
opposite sign. When pivoting from an ACL state s on that path to t — f{sj i) 
where l{si) is the duplicate label in r{s) = (^i, . . . ,Sm), then the CL state found 
at the end of the path by continuing to pivot in that direction has opposite sign to 
sign(5',/). There are as many CL states of sign 1 as of sign —1. 
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Proof. Assume that the pivoting system is oriented; otherwise complementary piv- 
oting (already described informally above) is part of the following description by 
disregarding all references to signs. Consider a CL state s and r{s) = {s\,. . . ^Sm), 
with w declared as the missing label for the path that starts at s, and let l{si) = w. 
We can define sign(5, /) as in (10), which is just sign(5) in (9), because l{r{s)) = 
(/(r(s)) I i^w). The following considerations apply in the same way if s is an ACL 
state with duplicate label l{si) . The path starts (or continues, if s is ACL) by pivoting 
tot = f{s, i) . Assume r^{t) = {r{s) | ? — )■ m) as in Definition 2. Then {l{r{s)) \ i ^w) 
is a permutation of [m] , which is equal to {l{r^{t)) | z — )> w) , and {l{r{t)) \ n{i) — > w) 
is a permutation of [m] with parity(;r) •parity(/(r(5)) \ i ^ w) as its parity. Hence, 
by (8) 

sign(s,/) = o{s) ■parity(/(r(5)) \ i ^ w) 

= —(y{t) ■parity(7r) ■ parity(/(r(5)) \i^w) 
= -(7{t) ■parity(/(r(r)) | 7c{i) w) 
= -sign{t,7t{i)). 

If l{u) is the missing label w, then t is the CL state at the other end of the path 
and sign(?) = sign(?, n{i)), which is indeed the opposite sign of the starting state s. 
Otherwise, label l{u) is duplicate, with l{u) = l{sj) for some j ^ i, that is, l{tji{{)) = 
K^n{j)) for r{t) = (?!,...,?„), so that the path continues with the next pivoting step 
from t to f{t, Tt{j)), where by (1 1) 

sign(?, 7t{j)) = -sign(?, 7J:{i)) = sign{s,i), 

that is, this step continues from a state with the same sign as the starting CL state, 
and the argument repeats. This proves the theorem. □ 

For a labeled poly tope P as in (1), an oriented pivoting system is obtained as 
follows: The states in S are the vertices x of P, and by the assumptions on P following 

(1) each vertex x lies on exactly m facets F^^ for p E [m], where we take r{x) = 
(^i, . . - jSm) as the representation of x with si,. . . ,Sf„ in any fixed order. Moreover, 
the normal vectors of these facets in (2) are linearly independent. For any i in 
[m] , the set flpe [«]-{;} ^sp is an edge of P with two vertices x and y as its endpoints, 
which defines the pivoting function as y = f{x, i) . The orientation of the vertex x is 
given by 

a{x) = sgn(det[a,i ■ --asj) (12) 

with the usual sign function sgn(z) for reals z and the determinant detA for any 
square matrix A. The following proposition is well known (Lemke and Grotzinger, 
1976, for example, argue with linear programming tableau entries; Eaves and Scarf, 
1976, Section 5, consider the index of mappings); we give a short geometric proof. 

Proposition 4 A labeled polytope P with orientation o{x) as in (12) for each vertex 
X of P defines an oriented pivoting system. 

Proof. Consider pivoting from x to vertex y = f{x, i) . We want to prove (8), that 
is, o{y) — —(y{x) •parity(;r) where n is the permutation so that r^{y) — {r{x) \ 
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/ — > m). Let X be on the m facets Fs^jFs2, ■ ■ ■ ,Fs^ as in (2). The representation 
r(x) = (i'l, . . .,5m) determines the order of the columns of the matrix [fljia^j ■ ■ -^.v^] 
whose determinant determines the orientation a{x) in (12). Any permutation of the 
columns of this matrix changes the sign of the determinant according to the parity of 
the permutation, so for proving (8) the actual order of (^i, . . .,5m) in r{x) does not 
matter as long as it is fixed. Hence, we can assume that 7t is the identity permuta- 
tion, and that pivoting affects the first column (/ = 1), so that y is on the m facets 

We show that det[(2vy(2v2 ■ ' ■'^.vml ^nd det[(2vi«,s2 • ' ■'^.s,,,] h^^^ opposite sign, that is, 
o{y) — —o{x) as claimed. The m + 1 vectors ajQ,aij,aj2 5 • ■■■>^sm are linearly depen- 
dent, so there are reals cq, ci , . . . , c^, not all zero, with 

m 

Lcp< = OT. (13) 

Note that cq ^ 0, because otherwise the normal vectors , a^j, . . . , «.s„ of the facets 
that define x would be linearly dependent, and similarly cj 7^ 0. Multiply the sum 
in (13) with both y and x, where aj^^y = aj^x = bg^ for p — 2, . . . ,m. This shows 

^0 <^jQy + c\ cil^ — ^0 clIq^ + ^1 ^s\^ equivalently 

co(«Jy - «J^) = ci (ajx - ajy) , 

so Co and c\ have the same sign because x is not on facet F^^ and y is not on facet 
F^j , so al^y - a^x = h^^ - aj^x > and aj^x - aj^y = bs^ - aj^y > 0. By (13), 

= det[{asoCo + ci ) • • • ^sj = cq det[asQas2 ■ ■ ■ agj + c\ detfa^iflij " " ' ^s,^ 

which shows that detfaiga^j ■ ■ -^im] and detfa^^a^j " ' "^^m] ^^^^ indeed opposite sign. 

□ 

The orientation of a vertex of a simple polytope P depends only on the determi- 
nant of the normal vectors aj of the facets in (12), but not on the right hand sides hj 
when P is given as in (1). Translating the polytope by adding a constant vector only 
changes these right hand sides. If is in the interior of P, then one can assume that 
hj = 1 for all j in [n] . The convex hull of the vectors aj is then a simplicial polytope 
P^ called the "polar" of P (see Ziegler, 1995). The vertices of P^ correspond to the 
facets of P and vice versa. A pivoting system for the simplicial polytope has its ver- 
tices as nodes and its facets as states, which one may see as a more natural definition. 
However, the facets of a simplicial polytope are oriented via (12) only if it has in its 
interior, which is not required for the simple polytope P. For common descriptions 
such as (4), we therefore prefer to look at simple polytopes. 

Theorem 3 and Proposition 4 replicate, in streamlined form, Shapley's (1974) 
proof that the equilibria at the ends of a Lemke-Howson path have opposite index. 
Applied to the polytope P in (4), the completely labeled vertex does not represent a 
Nash equilibrium, and it is customarily assumed to have index — 1 , which is achieved 
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by multiplying all signs with —1 if m is even. Lemke and Grotzinger (1976) have 
generalized Shapley (1974) by considering abstract manifolds, which are special piv- 
oting systems, proceeding very similarly to this section. We consider these manifolds 
as special cases of Euler complexes in the next section. 

3 Oriented Euler complexes 

Edmonds (2009) introduced the concept of an Euler complex or "oik" to abstract 
from path-following arguments such as the Lemke-Howson algorithm. Using the 
new concept of an oriented oik, this section shows how to give the path a direction, 
specifically for "room partitions" in an oik family. The room partitions normally 
have to be ordered. When all oiks in the family are a single oik, the unordered room 
partitions are known to be connected by paths. Then the paths can also be given a 
direction if the oik is oriented and has even dimension (Theorem 11). We discuss the 
connection of labels with "Spemer oiks" in Appendix A. 

Definition 5 Let V be a finite set of nodes and d be an integer, d >2. A J-dimen- 
sional Euler complex or d-oik on V is a multiset M of J -element subsets of V, called 
rooms, so that any set W of J — 1 nodes is contained in an even number of rooms. 

If W is always contained in zero or two rooms, then the oik is called a manifold. A 
wall is a (J — 1) -element subset of a room R. A neighboring room to R for W is any 
room that contains W. 

In the preceding definition we follow Edmonds, Gaubert, and Gurvich (2010) of 
choosing d rather than d—l (as in Edmonds, 2009) for the dimension of the oik. A 
2-oik on y is an Euler graph with node set V and edge multiset M. We allow for 
parallel edges (which is why M in Definition 5 is a multiset, not a set) but no loops. 

Rooms are often called "abstract simplices", and a longer term for manifold is 
"abstract simplicial pseudo-manifold" (e.g., Lemke and Grotzinger, 1976). The fol- 
lowing definition generalizes the common definition of coherently oriented rooms in 
manifolds (Hilton and Wylie, 1967, p. 54) to oiks. 

Definition 6 Consider a d-oik ^ on V and fix a linear order on V . Represent each 
room R = {si,.. ..sj} in M as r{R) = {si,...,sj) where i'l, . . . ,5^ are in increasing 
order. Choose an orientation o{R) in { — 1,1}. The induced orientation on any 
wall W = R — {si} is defined as (— 1)'<7(/?). The orientation of the rooms is called 
coherent, and the oik oriented, if half of the rooms containing any wall W induce 
orientation 1 on W and the other half orientation —1 on W. 

As an example, consider a 2-oik, where rooms are the edges of a Euler graph. 

Suppose an edge {m,v} is oriented so that a{u,v) = 1. Then the induced orientation 
on the wall {u} is —1 and on {v} it is 1, so {u,v} becomes the edge (m,v) of a 
digraph oriented from m to v. A coherent orientation means that each wall (that is, 
node) has as many incoming as outgoing edges, so this is an Eulerian orientation of 
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the graph (which always exists; for <i > 2 there are already manifolds that cannot be 
oriented, for example a triangulated Klein bottle). In general, the simplest oriented 
oik consists of just two rooms with equal node set but opposite orientation. In an 
oriented Euler graph, this is a pair of oppositely oriented parallel edges. 

Proposition 7 A d-oik M onV defines a pivoting system {S,V,m,r,f) as follows: 
Let S = m = d, and r and o be as in Definition 6. For any wall W , match 
the 2k rooms that contain W into k pairs {R,R'), where R and R' induce opposite 
orientation on W if the oik is oriented. Then f{R, i) — R' if r{R) — {si,.. .,Sd) and 
W = R — {si}. If G is coherent, then the pivoting system is oriented. 

Proof. Let RUR' — {^i, . . . — RU {sj} — R' U {si}, with ^i, . . . j^^z+i in in- 

creasing order, and let i < j, otherwise exchange R and R'. Then r{R') is obtained 
from r{R) by replacing Si with sj followed by the permutation 7t that inserts sj at 
its place in the ordered sequence by "jumping over" j — i—l elements . . . , sj-i 
to remove as many inversions, so parity(;r) = ( — 1)-'^'^^. Hence, f{R,i) = R' is 
well defined. If o is coherent, then R and R' induce on the common wall RDR' 
the opposite orientations (— l)'a(7?) and (— l)-'^a(i?') (because Si ^ R'), that is, 
a{R') = -a{R){-l)j-'-^ = - a {R) -parity {n) as requiredin (8). □ 

The matching of rooms with a common wall into k pairs described in Proposi- 
tion 7 is unique if the oik is a manifold. In a 2-oik, that is, an Euler graph, such a 

matching of incoming and outgoing edges of a node is for example obtained from an 
Eulerian tour of the graph, which also gives a coherent orientation. 

For an "oik-family" . . .,^h where each is a dp-oik on the same node 
set y for p e [h], Edmonds, Gaubert, and Gurvich (2010) define the "oik-sum" as 
follows. 

Definition 8 Let ^p be a dp-oik on V for p e [h] , and m — Lp=i dp. Then the oik- 
sum M = h is defined as the set of m -element subsets R of [h] x V so 

that 

R^Ri\SR2^---^Rh^{{l}xRi)U{{2}xR2)U---U{{h}xRh) (14) 

where Rp G Mp for p e[h]. For a fixed order < on V, we order [h] x V lexicograph- 
ically by {p, u) < {q, v) if and only if p < q, or p = q and u<v. 

As observed by Edmonds, Gaubert, and Gurvich (2010), the oik-sum ^ is an oik. 
A neighboring room of /? = i?i l+li?2 W ■ ■ ■ W^/i is obtained by replacing, for some p, 
the room Rp with a neighboring room R'p in Mp . The next proposition states, as a new 
result, that the oik- sum is oriented if each ^p is oriented. According to Definition 6, 
this requires an order on the node set [h] x V to yield an order on the nodes in room 
R in (14), which is provided in Definition 8: The nodes of each room Rp are listed 
in increasing order (on V), and these J^-tuples are then listed in the order of the 
rooms Ri,...,Rii; this becomes the representation r{R) used to define the orientation 
o on ^. 
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Proposition 9 The oik-sum ^ in Definition 8 is an m-oik over [h\ x V . If each 
is oriented with Op, so is 3$, with 

h 

p=\ 

Proof. Clearly, each room of =^ as in (14) has m elements. Any wall of i? is 
given by W = — { (p, v) } for some p in [h] and v in i?^ . Then any neighboring room 
R' in 3^ for W is given by 

R' = Rx^- ■\SRp-\ l±li?p • • • l±l/?/, 

for the neighboring rooms R'p in Mp for Rp — {v}, of which, including Rp, there is 
an even number. This shows that M is an m-oik. 

For the orientation of M if each Mp is oriented with Op, represent R as r{R) by 
listing the elements of R in lexicographic order as in Definition 8. Then the induced 
orientation on any wall W = R — {{p,v)} as in Definition 6 is obtained from the 
induced orientation on Rp — {v}, as follows. Suppose 5p . . . ,5^ are the nodes in 

Rp in increasing order, where v = . Then the induced orientation on Rp — {v} in 
Mp is {—iyOp{Rp). In r(i?), node v appears in position Lj=i dj + / , so the induced 
orientation of on W is, with (y{R) is defined as in (15), 

(-I)^^^1'^^+'CT(^) = (-I)'c7p(i?p)(-1)^&'^^- n ^'M- (16) 

qe[h]-p 

All the rooms in ^ that contain W are obtained by replacing Rp with any room 
R'p that contains Rp — {v}. Half of these have induce the same orientation as Rp 
on Rp — {v}, half of these the other orientation. Because this affects only the term 
{ — \yOp{Rp) in (16), half of the rooms R' that contain W induce one orientation on 
W and half the other orientation. So a is a coherent orientation of □ 

Consider now an oik-family .. ,^%h where S?,p is a J^-oik on V for p in [h] 
so that |y I = m~ dp. Suppose Rp e ^p for p in [h] and [jp^i Rp — V (so the 
rooms Rp are, as subsets of V, also pairwise disjoint). Then (^i, . . . ,i?/,) is called an 
ordered room partition. In the following theorem, the even number of ordered room 
partitions is due to Edmonds, Gaubert, and Gurvich (2010); the observation on signs 
is new. 

Theorem 10 Let Mp be a dp-oik on V for p in [h] and \ V\ = m = Lp=i dp. Then the 
number of ordered room partitions is even. If each ^p is oriented as in Proposition 9, 
then there is an equal number of ordered room partitions of positive as of negative 
sign, where the sign of a room partition {Ri,. . . ,Rh) is defined by 

sign(i?) = sign(i?i, . ..,Rh)^ <y{Ri W . . . Wi?^) •parity(;r) (17) 

with the permutation % ofV given according to the order of the nodes ofV in r{R), 
that is, with 7t{u) < n{v) if u & Rp and v & Rq and p < q, or u,v & Rp and u <v 
inV. 
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Proof. This is a corollary of Theorem 3 and Propositions 7 and 9. Assume that 

V — {vi,.. . , Vot} with the order on V given by v,- < vj for z < j (or just let V = [m]). 
Define the labeling / : [/z] x V — )> [m] by l{p,Vi) = / for / G [m] . Then the CL rooms 

i?i l±l . . . I±l7?/, of ^1 H h=^/j are exactly the ordered room partitions, with the sign 

in (17) defined as in (9). So there is an equal number of them of either sign. 

If the oiks are not all oriented, then the paths that connect any two CL states are 
still defined, so the number of ordered room partitions is even, except that they have 
no well-defined sign. □ 

Connecting any two room partitions by paths of ACL states as in the preceding 
proof corresponds to the "exchange graph" argument of Edmonds (2009), where the 
ACL states correspond to skew room partitions . . . ,i?/,) defined by the property 
Up=i Rp — V — {w} for some w in V ; here w represents the missing label. 

Suppose now that all oiks in the oik family are the same J-oik M' over V 
for p in [h] , with \V\ = m = h - d. Then any ordered room partition , . . . , i?/,) de- 
fines an (unordered) room partition . . . ,Rh}- Any such partition gives rise to 
h\ ordered room partitions, so if > 2 their number is trivially even. However, the 
path-following argument can be applied to the unordered partitions as well (which is 
the original exchange algorithm of Edmonds, 2009), which shows that the ordered 
room partitions at the two ends of the pivoting path define different unordered parti- 
tions. The next theorem shows that unordered partitions {^i, . . . are connected 
by pivoting paths, which are essentially the same paths as in Theorem 10, and that 
the sign property continues to hold when d is even and is oriented. 

Theorem 11 Let he a d-oik on V and \V\ = m = h ■ d. Then the number of 
room partitions . . . is even. If is oriented with o' and d is even, then 
sign{Ri, ... ,R}i) as defined in (15) with Op — o' and (17) is independent of the 
order of the rooms Ri,. . . ,Rfi, and there are as many room partitions of sign 1 as of 
sign —1. 

Proof. We consider unordered multisets {R\_,. . . ,Rh} of h rooms of as states s 
of a pivoting system. We first define a representation r{s) = (51, . . .,5m)- Let Rp = 
{i'j', . . . , 5^} for /7 in [h] where if , . . . , 5^ are in increasing order according to the order 
on V . Fix some order of the rooms in M' , for example the lexicographic order with 
some tie-breaking for rooms that have the same node set. Assume that the rooms 
R\,...,Rh are in ascending order, which defines a unique representation of s as 

r{s) = r{{R,,...,Rh}) = {si...,s]i,sl...,s% . . . ,5?, . . . ,4). (18) 

(Note that r may not be injective, which is allowed.) Assume that neighboring rooms 
in ^' are matched into pairs Rp^R'p containing the wall Rp — {si\ as in Proposition 7. 

The pivoting step from s = {R\, . . .^Rh} to t = f{s, i) replaces Rp by R'^. 

In (18), the nodes of each individual room Rp still appear consecutively as in 
the permutation n in Theorem 10, except for the order of the rooms themselves. 
Then with vi , . . . , as the nodes of V in increasing order and the "identity" labeling 
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/ : y — )■ [m], l{vi) = i, the m-tuple l{r{s)) defines a peraiutation n of [m] if 5 is a 
room partition, as in (17). Then the parity of n does not depend on the order of the 
rooms in 5 if J is even, so the sign in (17) is well defined and the same as in (9). An 
ACL state 5 is a skew room partition, which has two opposite signs as in (1 1). Then 
the claim follows from Theorem 3. □ 




Figure 1: A 3-oik with triangles as rooms. The circular arrows indicate the positive 
orientation of nodes in a room. 

The following example shows that we cannot expect to define a sign to unordered 
room partitions when M' has odd dimension d (see also Merschen, 2012, Figure 3.6). 
Let d = 3 and consider the oik defined by the eight vertices of the 3 -dimensional 
cube, which correspond to the facets of the octahedron, shown as the triangles in Fig- 
ure 1 including the outer triangle marked "A". A coherent orientation of the eight 
rooms is obtained as follows (shown in Figure 1 with a circular arrow that shows the 
positively oriented order of the nodes): o{A) = a{l23) = 1, o{B) = a{l45) = — 1, 
(t(C) = a(124) = -1, o{D) = a(135) = 1, G{a) = C7(456) = 1, o{b) = (7(236) = 
— l,o{c) — a (356) = — 1, o{d) — a (246) = 1. The four room partitions are {A, a}, 
{B, b}, {C, c}, {D, d} . Any two of these are connected by pivoting paths, so they can- 
not always have opposite signs at the end of these paths. However, for ordered room 
partitions the signs work. For example, (A, a) is connected to {b,B) via the com- 
plementary pivoting steps (123,456) —?■ (236,456) — > (236, 145), and to (C,c) via 
the steps (123,456) (124,456) (124,356). Moreover, (C,c) connects to {B,b) 
via (124,356) ^ (145,356) ^ (145,236). We have sign(A,a) = 1, sign{b,B) - -1 
(because 236145 has parity —1), and sign(C,c) = — 1 and sign(5,&) = 1. The two 
ordered room partitions {b,B) and (5, b) have different signs because they define two 
permutations 236 145 and 145236 of opposite parity. 

If the oik cannot be oriented, then Lemke and Grotzinger (1976) have shown 
(for nonorientable manifolds) that opposite signs for CL rooms cannot be defined in 
general; see also Grigni (2001). 

We discuss the close connection of labels and "Sperner oiks" considered by Ed- 
monds, Gaubert, and Gurvich (2010) in Appendix A. There we argue that the defi- 
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nition of "sign" requires a reference to the permutation of the labels which does not 
seem simpler when looking at room partitions with a Spemer oik instead. 

4 Signed perfect matchings 

This section is concerned with algorithmic questions of room partitions in 2-oiks, 
which are perfect matchings in Euler graphs. The sign of a perfect matching, for 
any orientation of the edges of a graph, is closely related to the concept of a Pfaffian 
orientation of a graph, where all perfect matchings have the same sign. The com- 
putational complexity finding such an orientation is an open problem (see Thomas, 
2006, for a survey). An Eulerian orientation is not Pfaffian by Theorem 1 1, a fact that 
is also easy to verify directly. The main result of this section (Theorem 12) states that 
in a graph with an Eulerian orientation, a second perfect matching of opposite sign 
can be found in polynomial (in fact, near-linear) time. This holds in contrast to the 
complementary pivoting algorithm, which can take exponential time; Casetti, Mer- 
schen and von Stengel (2010) have shown how to apply results of Morris (1994) for 
this purpose. However, the pivoting algorithm takes linear time in a bipartite Euler 
graph, and a variant can be used to find an oppositely signed matching in a bipartite 
graph that has no source or sink (Proposition 13). 

We follow the exposition of Pfaffians in Lovasz and Plummer (1986, Chapter 8). 
The determinant of an m x m matrix B with entries bij is defined as 

m 

det5 = J^parity(;r)n^;r(0 (19) 

7t i=l 

where the sum is taken over all permutations k of [m] . Let B be skew symmetric, 
that is, B = -B^. Then det5 = det{-B'^) = det(-5) = (-l)"det5, so det5 = if 
m is odd. Assume m is even. Then 

det5=(pf5)2 (20) 

for a function pf B called the Pfaffian of B, defined as follows. Let ^{m) be the set 
of all partitions s of [m] into pairs, s = {{51,52}, . . . , {sm-iiSm}}, and let parity(5) be 
the parity of (51,52, • • • ,>5'm) seen as a permutation of [m] under the assumption that 
each pair {52yt-i , S2k} is written in increasing order, that is, 52^-1 < 52^ for kin [m/2] ; 
the order of the pairs themselves does not matter. Then 

m/2 

pf5= parity 

(21) 

s€^{m) k=l 

In fact, because B is skew symmetric, the order of a pair {s2k-i,S2k) can also be 
changed because this also changes the parity of 5. An example of (21) is m — 4 
where pf 5 = &12&34 — ^13^24 + ^14^23 • 
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Parameswaran (1954) and Lax (2007, Appendix 2) show that a skew-symmetric 
matrix B fulfills (20) for some function pf 5. For a direct combinatorial proof, one 
can see that the products in (19) are zero for those permutations k where n{k) =k for 
some k, and cancel out for the permutations with odd cycles; then only permutations 
with even-length cycles remain, which can be obtained uniquely via pairs of partitions 
taken from ^{m) (see also Jacobi, 1827, pp. 354ff, and Cayley, 1849). 

Consider a simple graph G with node set [m] . An orientation of G creates a 
digraph by giving each edge {m, v} an orientation as (m, v) or (v, u) . Define the m x m 
matrix B via 

{0 if {m, v} is not an edge, 
1 if {m, v} is oriented as (m, v) , (22) 
— 1 if {m, v} is oriented as (v, u). 

Then B is skew symmetric. Any s in ^{m) is a perfect matching of G if and only if 

Y[k=i^s2k-i,s2k 7^ 0' ^° ^'^^y perfect matchings of G contribute to the sum in (21). 

If G is an Euler graph, that is, a 2-oik with edges as rooms, and the orientation is 
Eulerian so that every node has equal in- and outdegree, then this defines the orien- 
tation of edge {m, v}, assuming m < v, as (7({m,v}) = buv, according to Definition 6. 
Then by (15) and (17), a perfect matching s has the sign 

m/2 

sign(s) = parity(5i, . . . , ■ ^^2^-1,^2* ' 

so the Pfaffian pf 5 in (21) is the sum over all matchings of G weighted with their 
signs. For the Eulerian orientation, that sum is zero by Theorem 11, which follows 
also from (20) because Bl = 0, so det B = 0. 

In our Definition 5 of a J-oik, M can be a multiset, which for d = 2 defines 
an Euler graph G which may have parallel edges and then is not simple. The rooms 
themselves have to be sets, so loops are not allowed. In this case, (22) can be extended 
to define buv as the number of edges oriented as (m,v) minus the number of edges 
oriented as (v, u) . This counts the number of matchings with their signs correctly; 
oppositely oriented parallel edges (m, v) and (v, u) cancel out both in contributing to 
buv and when counting matchings with their signs. 

For any graph G and any orientation of G, the sign of a perfect matching s is 
most easily defined by writing down the nodes of each edge {s2k-uS2k} in the way 
the edge is oriented as (52/t-i , S2k) ; this does not affect (21) as remarked there. When 
writing down the nodes si,...,Sm this way, sign(5) = parity (51 , . . . , Sm) and pfB — 
Y^seJ^iG) sign(5) where M{G) is the set of perfect matchings of G. 

A Pfaffian orientation is an orientation of G so that all perfect matchings have 
positive sign. Its great computational advantage is that it allows to compute the num- 
ber of perfect matchings of G using (20) by evaluating the determinant det 5, which 
can be done in polynomial time. In general, counting the number of perfect matchings 
is #P-hard already for bipartite graphs; the question if a graph has a Pfaffian orienta- 
tion is polynomial-time equivalent to deciding whether a given orientation is Pfaffian 
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(see Vazirani and Yannakakis, 1989, and Thomas, 2006). For bipartite graphs, this 
problem is equivalent to finding an even-length cycle in a digraph, which was long 
open and shown to be polynomial by Robertson, Seymour, and Thomas (1999). For 
general graphs, its complexity is still open. 

We now consider the following algorithmic problem: Given a graph with an Eule- 
rian orientation and a perfect matching, find another matching of opposite sign, which 
exists. Without the sign property, a second matching can be found by removing one 
of the given matched edges from the graph and applying the "blossom" algorithm of 
Edmonds (1965) to find a maximum matching, which finds another perfect matching 
for at least one removed edge; however, its sign cannot be predicted, and adapting 
this method to account for the sign seems to lead to the difficulties related to Pfaffian 
orientations in general graphs. Merschen (2012, Theorem 5.3) has shown how to find 
in polynomial time an oppositely signed matching in a planar Euler graph, and his 
method can be adapted to graphs that, like planar graphs, are known to have a Pfaffian 
orientation. 

The following theorem presents a surprisingly simple algorithm for any Euler 
graph. It runs in near-linear time in the number of edges of the graph and is faster 

and simpler than using blossoms. The inverse Ackermann function a is an extremely 
slowly growing function with a(n) < 4 for n < 2^^^^ (Gormen et al., 2001, Section 
21.4). 

Theorem 12 Let G = {V,E) be a graph without loops with an Eulerian orientation, 
and let M he a perfect matching of G. Then a perfect matching M' of opposite sign 
can be found in time 0{\E\- a{\V\)), where CC is the inverse Ackermann function. 

Proof. The matching M is a subset of E. A sign-switching cycle C is an even-length 
cycle (as a set of edges) so that exactly every other edge in C belongs to M and 
so that C has an even number of forward-oriented edges (in either direction of the 
cycle). Then the symmetric difference M' = MAC has opposite sign to M. To see 
this, suppose first that all edges in C point forward, and that COM consists of the first 
k/2 edges {si,S2}, ■ ■ ■ , (sk-uSk) of M (which does not affect the sign of M). Then 
these edges are replaced in M' by (sk.si), (52, ■^s), • • • , (sk-ii^k-i), which defines 
an odd permutation of these k nodes, so M' has opposite sign to M. Ghanging the 
orientation of any two edges in C leaves the sign of both M and M' unchanged (if 
both edges belong to M or to M') or changes the signs of both M and M', so they 
stay opposite. This proves the claim. 

So it suffices to find a sign- switching cycle C for M, which is achieved by the 
following algorithm: Successively apply one of the following reductions (a) or (b) to 
G until (c) applies: 

(a) If V in y has indegree and outdegree 1 with edges (w, v) and (v, w) , then if w = w 
go to (c), otherwise remove v from V and (m,v) and (v,w) from E and contract 
u and w into a single node. 

(b) If Z) is a directed cycle of unmatched edges (so D cE —M), remove all edges in 
D from E. 
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(c) The two edges (m, v) and (v, u) , one of which is matched, form a sign-switching 
cycle C of the reduced graph. Repeatedly re-insert the edge pairs (w', v') , (v', w') 
removed in the contraction (a) into C until C is a cycle of the original graph. 
Return C. 

Steps (a) and (b) preserve the invariant that G has an Eulerian orientation and a perfect 
matching. Namely, in (a) one node and one matched and one unmatched edge is 
removed from G, and the two contracted nodes u and w together have the same in- 
and outdegree and an incident matched edge. In (b), all nodes of the cycle D have 
their in- and outdegree reduced by 1 . If reduction (a) cannot be applied because every 
node has at least two outgoing edges, then one of them is unmatched, and following 
these edges will find a cycle D as in (b). So the reduction steps eventually terminate. 
In each iteration in (c), the two re-inserted edges {u',v') and (v',w') point in the 
same direction and one of them is matched, so this preserves the property that C is 
sign-switching. 

The above algorithm is clearly polynomial. The Appendix describes a detailed 
implementation with near-linear running time in the number of edges, and gives an 
example. Its essential features are the following. The algorithm starts with the end- 
point of a matched edge, and follows, in forward direction, unmatched edges when- 
ever possible. It thereby generates a path of nodes connected by unmatched edges. 
If a node is found that is already on the path, then some final part of that path forms 
a cycle D of unmatched edges that are all discarded as in (b). Then the search starts 
over from the beginning of the cycle that has just been deleted. If, in the course of 
this search, a node v is found where the only outgoing edge (v, w) is matched, then 
the contraction in (a) applies with (m, v) as unmatched edge. The matched edge (v, w) 
is remembered as the original matched edge incident to w, with (m, v) as its "partner", 
for possible later re-use in (c). The two edges are removed from the lists of incident 
edges to u and w. Edges are stored in doubly-linked lists that can be moved and 
deleted from in constant time. The endpoint w of the matched edge {u,w) contracted 
in step (a) may be a node that has been visited on the path, so that the reduction (b) 
immediately follows; if w is the first node of the path, the search has to re-start. 

Contracted nodes of the reduced graph are represented by equivalence classes of 
a standard union-find data structure, which can be implemented with amortized cost 
a(|y I) per access (Tarjan, 1975). Contracting u and w in (a) is done by applying the 
"union" operation to the equivalence classes for u and w, and any node is represented 
via the "find" operation applied to an original node. The nodes in edge lists are 
always the original nodes, so that each edge is visited only a constant number of 
times, resulting in the running time 0(|£'|-a(|y|)). 

As described in Appendix B in Figure 10, the cycle C in (c) is obtained by recur- 
sively re-inserting matched edges (v',w') and their "partners" {u',v") until the nodes 
v' and v" do not just belong to the same equivalence class (as at the time of contrac- 
tion) but are actually the same original node, v' = v", of G; a similar recursion is 
applied to the other nodes u' and w'. □ 
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In the remainder of this section, we consider the complementary pivoting algo- 
rithm for perfect matchings in Euler graphs. Given a graph G with an Eulerian ori- 
entation and a matching M, the complementary pivoting algorithm finds a matching 
with opposite sign to M. The intermediate steps of the algorithm use skew match- 
ings, which are sets of edges that cover all nodes except for one missing node w, and 
that are pairwise disjoint except for one pair of edges that intersect in one duplicate 
node u. This is a special case of the pivoting system used to prove Theorem 11, with 
matchings as CL states and skew matchings as ACL states with missing label w and 
duplicate label u. In "oik terminology", the rooms are here edges, the walls nodes, 
and (skew) room partitions are (skew) matchings. 

The algorithm uses a pairing of the k incoming edges to the k outgoing edges of 
a node as in Proposition 7 to pivot from one edge to the next. Such a pairing defines 
a set of Eulerian cycles that cover all edges (which may be a single Euler tour of the 
graph), and vice versa. In a pivoting step, a matched edge {m,v} is replaced by its 
paired edge {v,m'}, where these edges point in the same direction either as (m,v), 
(v,m') or as (v,m), (m',v). When starting from a perfect matching, u is the missing 
node w. In a skew matching, u is the duplicate node. The pivoting step reaches either 
a new skew matching where u' is the duplicate node, or a perfect matching (which 
terminates the path) when u' = w. 

If G is bipartite, then this algorithm terminates in time 0(| V|), as noted by Mer- 
schen (2012, Lemma 4.3). In fact, a simple extension of the pivoting method applies 
to general bipartite graphs which are oriented so that the graph has no sources or 
sinks (which shows that such an orientation is not Pfaffian). 

Proposition 13 Consider a bipartite graph G = {V,E) with an orientation so that 
each node has at least one incoming and outgoing edge, and a perfect matching M. 
Then a matching of opposite sign can be found in time 0{\V\). 

Proof. The algorithm computes a path of nodes mo,mi, . . . until that path hits itself 
and forms a cycle C, which will be sign-switching with respect to M. The edges 
on the path are successive matched-unmatched pairs of edges {"i/ti M2/t+i} in M 
and {u2k+iTU2k+2} h\ E — M iox k>0 that point in the same direction either as 
{u2k:U2k+i), {u2k+i:U2k+2) or as (m2A:+i , M2A:) , ( M2;t+2 , M2yk+ 1 ) • Starting from any nodc 
Mo and k = 0, these are found by following from node U2k its incident matched edge 
to U2k+i, where this node has an outgoing unmatched edge to U2k+2 in the same di- 
rection because U2k+i has at least one incoming and one outgoing edge. This repeats 
with k incremented by one until U2k+2 is a previously encountered node, which is 
of the form U2i for some < i < k because the graph is bipartite. Then the nodes 
M2i, • • • , U2k+2 define a cycle C which is sign-switching because it has an even number 
of forward-pointing edges. Hence, MAC is a matching of opposite sign to M. Each 
node is visited at most once, so the running time is 0(|y |) . □ 

If G is not bipartite, then the complementary pivoting algorithm may have ex- 
ponential running time, for any starting node that serves as a missing label. The 
construction is adapted from the exponentially long Lemke paths of Morris (1994) 



20 



for labeled dual cyclic polytopes. The completely labeled vertices of such polytopes 
correspond to perfect matchings in Euler graphs, as noted by Casetti, Merschen, and 
von Stengel (2010), in the following way. 

A dual cyclic polytope is defined in any dimension m with any number n of facets, 
n > m, as the "polar polytope" of the convex hull of n points }i{tj) on the moment 
curve = (?,?^, ...,?"')^ for j in [n] (see Ziegler, 1995). Its vertices have been 
described by Gale (1963): The m facets that a vertex x lies can be described by a bit 
string g = gig2 ■ ■ - in {0, 1}" so that gj = 1 if and only if x is on the jth facet, for 
j in [n] . Then these bit strings fulfill the evenness condition that whenever g has a 
substring of the form 01^0, then k is even. We consider even m, so that these strings 
are preserved under cyclical shifts. The set G(m, n) of these "Gale strings" encodes 
the vertices of the polytope, and pivoting, and an orientation, can be defined in a 
simple combinatorial way on the strings alone. 

With a labeling /:[«]—)■ [m], the CL Gale strings therefore come in pairs of 
opposite sign. They correspond, including signs, to the perfect matchings of the graph 
G with node set [m] and (oriented) edges (/(j), /(y + 1)) for 1 < 7 < n and (/(n), /(I)) 
(Casetti, Merschen, and von Stengel, 2010; Merschen, 2012, Theorem 3.4). That is, 

the cyclic sequence /(I), /(n), /(I) defines an Eulerian orientation of G, and G is 

an Euler graph. The graph has parallel edges and possibly loops, where the latter can 
be omitted. The I's in a Gale string come in pairs, which correspond to edges of G. A 
pivoting step from one ACL Gale string to another means that a substring of the form 
1^*^0 is replaced by 01^*^, which translates to k pivoting steps of skew matchings in G. 
Morris (1994) gives a specific labeling for n = 2m where all complementary pivoting 
paths, for any dropped label, are exponentially long in m. The corresponding Euler 
graph and the pivoting steps are described in Merschen (2012, Section 4.4). 

5 Conclusions 

We conclude with open questions on the computational complexity of pivoting sys- 
tems. 

Consider a labeled oriented pivoting system whose components (in particular the 
pivoting operation) are specified as polynomial-time computable functions. Assume 
one CL state is given. Then problem of finding a second CL state belongs to the 
complexity class PPAD (Papadimitriou, 1994). This problem is also PPAD-complete, 
because finding a Nash equilibrium of a bimatrix game is PPAD-complete (Chen and 
Deng, 2006), which is a special case of an oriented pivoting system by Proposition 1 . 
However, there should be a much simpler proof of this fact because pivoting sys- 
tems are already rather general, so that it should be possible to encode an instance 
of the PPAD-complete problem "End of the Line" (see Daskalakis, Goldberg, and 
Papadimitriou, 2009) directly into a pivoting system. 

Finding a Nash equilibrium of a bimatrix game is PPAD-complete, and Lemke- 
Howson paths may be exponentially long. Savani and von Stengel (2006) showed 
this with games defined by dual cyclic polytopes for the payoff matrices of both 
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players, and a simpler way to do this is to use the Lemke paths by Morris (1994). 
One motivation for the study of Casetti, Merschen, and von Stengel (2010) was the 
question if finding a second completely labeled Gale string is PPAD-complete. This 
is unlikely because this problem can be solved in polynomial time with a matching 
algorithm. For the complexity class PPADS, where one looks for a second CL state 
of opposite sign (Daskalakis, Goldberg, and Papadimitriou, 2009), this problem is 
also solvable in polynomial time with our algorithm of Theorem 12. 

However, for room partitions of 3-oiks, already manifolds, finding a second room 
partition is likely to be more complicated. Is this problem PPAD-complete? We leave 
these questions for further research. 
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Appendix A: Labeling functions and Sperner Oiks 

One of the original motivations to consider room partitions for oiks , . . . , with 
possibly different dimensions is to abstract from the original Lemke-Howson algo- 
rithm for possibly non-square bimatrix games, which alternates between two poly- 
topes, represented by ^\ and ^2 (Edmonds, 2009). Similarly, our proof of Theo- 
rem 3 shows complementary pivoting as an alternating use of the pivoting function 
and the labeling function. Edmonds, Gaubert, and Gurvich (2010) cast the use of 
labels (or "colors") in terms of room partitions with a special manifold called a 
Sperner oik. If / : V ^ [m] is a labeling function, then the rooms of the Sperner oik 
are the complements of completely labeled sets, that is, 

^o = {ecy I |2| = |\/|-m, l{y-Q) = \m\}. (23) 

This is a manifold because is a wall of a room Q of ^0 if and only if V — W has 
m-\-\ elements of which exactly two have the same label, so adding either element 
to W defines the two rooms that contain W . In addition to ,^o> suppose that M is an 
m-oik on Y and defines a pivoting system as in Proposition 7. Then an ordered room 
partition (i?, Q) with R e ^ and Q E Mq is just a completely labeled room R of M. 
Complementary pivoting with missing label w amounts to the "exchange algorithm" 
with skew room partitions, which are our ACL states. 

Is the use of room partitions where one room comes from a Sperner oik more 
natural than the concept of completely labeled rooms? Obviously, the definitions are 
nearly identical, but apart from that we want to make two comments in favor of using 
labels. 

First, Edmonds, Gaubert, and Gurvich (2010) note that a Sperner oik is "poly- 
topal", that is, its rooms correspond to the vertices of a simple polytope. They leave 
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the construction of such a polytope as an exercise, which we give here to show the 
connection to the unit- vector games in Proposition 1. 

Proposition 14 Let \V\ — {vi, . . . , Vn} and I :V [m] so that l{vi) = ifor i e [m]. 
Consider the m X {n — m) matrixA= [^/(v„+i) " ' '^/(vn)] withA^ = [a\---arn\ and 

Po = e M"-'" \Ay<\,y> 0}. (24) 

Then Pq is a simple polytope, and y is a vertex of Pq if and only if it lies on n — m 
facets and the m non-tight inequalities in (24) fulfill 

{i e [m] I a^^ < 1 } U {l{vm^j) I yj > 0} = [m] . (25) 

Proof. For each / in [m] let 

L{i) = e [r - m] I l{vm+j) = «}• (26) 

Then the /th row of Ay < 1 says ajy = L;eL(i)yy ^ 1- Let y G Pq- For each /, if 
ajy = 'LjeL{i)yj = 1' thenjy > for at least one j in L(i), so i e {l{Vfn+j) \ yj > 0}, 
which shows (25). 

The non-empty sets L{i) form a partition of [n — m], and if L{i) is empty then 
flj = and the inequality ajy < 1 is redundant. Therefore the inequalities (24) can 
be re-written as 

£ yj<h yj>0 UeL{i)). (27) 

For each i in [m] , (27) defines a simplex whose vertices are the unit vectors and in 
]rI^(OI (if /^(/) is empty, this is the one-point simplex {()}). Hence, Pq is the product 
of these simplices and therefore a simple polytope, so any vertex y of Pq is on exactly 
n — m facets. □ 

Proposition 14 can be applied to any Spemer oik of dimension n — m obtained 

from / : y — )■ [m] which has at least one room, taken to be {v^+i, . . . , v,,} by num- 
bering V suitably. The n inequalities in (24) have labels 1, . . . ,m,/(vm-|-i), . . . ,/(vot) ; 
they define facets of Pq except for redundant inequalities ajy < 1 where a, =0. 
Then the n — m tight inequalities for each vertex y of Pq define a room of because 
the labels for the m non-tight inequalities for y are the set [m] according to (25), in 
agreement with (23). 

Suppose is an m-oik given by the vertices of the polytope P in (4), with labels 
l,...,m,/(m+ l),...,/(n) for its n inequalities (the same labels as for Pq). Then 
an ordered room partition i?, Q with R e M and 2 G .^o is a completely labeled 
room R, or vertex x of P, with Q corresponding to a vertex y of Pq. Except for the 
vertex pair (0, 0) , this is a Nash equilibrium {x, y) of the unit- vector game (A, C^) in 
Proposition 1 . In that game, there is no reference to labels, which are encoded in the 
payoff matrix A that defines Pq, just as the labels are encoded in the rooms of =^o- 
Like unit vector games, Spemer oiks may offer a useful perspective, but we do not 
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Figure 2: Two oriented Euler graphs which show that the parity of the permutation of 
all nodes matters. 

think it is deep; moreover, they only have a simple structure as products of simplices 
described in (27). 

Secondly, Sperner oiks are oriented, and the labels used in the proof of Theorem 
10 and 1 1 are simply the nodes of V . Perhaps using a Sperner oik, rather than labels, 
may avoid referring to the parity of l{r{s)) for a room partition s as in (9) when 
defining the sign of 5? The following example shows that already when 5 is a room 
partition for a 2-oik, one has to refer to the parity of l{r{s)) in some way. Figure 2 
shows two cases of 2-oiks over V = {1,2,3,4} with an orientation. The left 
oik has the two room partitions {12,34} and {14,23}, where C7i(12) = (72(34) = 1 
and ai(14) = -1, a2(23) = 1. According to (15), this implies a(12,34) = 1 and 
(7(14,23) = —1, so the two room partitions have opposite orientation (it suffices to 
consider unordered room partitions because d is even, as noted in Theorem 11). 

Similarly, the right oik in Figure 2 has the two room partitions {12,34} and 
{13,24}, where C7i(12) = C72(34) = 1 and ai(13) = C72(24) = 1, so all orientations 
are positive and a(12,34) — 1 and a(13,24) = 1, so these two room partitions have 
equal orientation. The difference is that the room partition 5 = {14, 23} defines an 
even permutation l{r{s)) = (1,4,2,3) of V, whereas {13,24} defines the odd permu- 
tation (1,3,2,4). So the sign of a room partition has to refer to the order in which the 
labels appear. 

We think that labeled pivoting systems are a general and useful way of represent- 
ing path-following and parity arguments, certainly for complementary pivoting and 
room partitions in oiks. 

Appendix B: Implementation Details of Finding a Sign- 
Switching Cycle in an Euler Graph 

Theorem 12 states that an oppositely signed matching in a graph with an Eulerian 
orientation can be found in near-linear time in the number of edges. In this appendix, 
we describe the details of the implementation of the algorithm outlined in the proof 
of Theorem 12. 

When e is an edge from u to v, then we call u the tail and v the head of e, and 
both u and v are called endpoints of e. 
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The algorithm appUes reductions (a) and (b) to the graph until it has a trivial sign- 
switching cycle which is expanded as in (c) to form a sign-switching cycle of the 
original graph. The algorithm starts with a node that is the head of a matched edge, 
and follows, in forward direction, unmatched edges whenever possible. It thereby 
generates a path of nodes connected by unmatched edges. If a node is found that is 
already on the path, then some final part of that path forms a cycle D of unmatched 
edges that are all discarded as in (b). Then the search starts over from the beginning 
of the cycle that has just been deleted. 

If, in the course of this search, a node v is found with the only outgoing edge 
being matched, the contraction in (a) is performed as follows. Suppose the three 
nodes in question are m,v,w with unmatched edge e from m to v and matched edge 
m from v to w, and no other edge incident to v. We take the edges e and m and 
node V out of the graph and contract the nodes u and w into a single node (with 
the method SHRlNK(e,m) discussed below), which creates a reduced version of the 
graph. Throughout the computation, the current reduced graph is represented by a 
partition of the nodes with a standard union-find data structure (Tarjan, 1975). We 
denote by [x] the partition class that contains node x, which has as its representative a 
special node called FIND(jc), where find is one of the standard union-find methods; 
we usually denote a representative node with a capital letter. That is, any two nodes x 
and y are equivalent (in the same equivalence class) if and only if find(x) = find(};) . 
In the reduced graph, every edge is only incident to the representative FIND(x) of a 
partition class, and the information for nodes that are not representatives is irrelevant. 
Initially, all partition classes are singletons {jc}, which is achieved by calling the 
makeset(.^) method. The method UNITE (jc,};) for nodes x,y merges [jc] and \y\ into 
a single set. 

Figure 3 shows an implementation of these methods as in Cormen et al. (2001, 
Section 21.3). (In this pseudo-code, an assignment such as Z,F •<— x,y assigns x 
to X and y to F, so for example x, j would exchange the current values of x 

and y.) Each partition class is a tree with x.parent pointing to the tree predecessor 
of node x, which is equal to ;c if jc is the root. For this root, x.rank stores an upper 
bound on the height of the tree. The unite method returns the pair X, F of former 
representatives of the two partition classes, where X is the new representative of 
the merged partition class and F is the representative no longer in use, which we 
need in order to move edge lists in the graph. With the "rank heuristic" used in the 
UNITE operation and the "path compression" of the recursive find method, the trees 
representing the partitions are extremely flat, with an amortized cost for the find 
method given by the inverse Ackermann function that is constant for all conceivable 
purposes (see Tarjan, 1975, and Cormen et al., 2001, Section 21.3). 

Every node x of the graph has its incident edges stored in an adjacency list, which 
for convenience is given by separate lists x.outlist and x.inlist for unmatched outgoing 
and incoming edges, respectively, and the unique matched edge x.matched which is 
either incoming or outgoing. Every edge e is stored in a single object that contains the 
following links to edges: e.nextout, e.prevout, e.nextin, e.previn, which link to the 
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makeset(jc): 

x.parent x 
x.rank ■<— 

UNITE 

Z , y FIND (x) , FIND iy) 

if x.rank > Y.rank then 
Y. parent •<— X 
return X,F 

else 

X.parent <(— 7 
if X .rank = Y .rank then 
Y.rank ■<— Y.rank+ 1 
return F,X 

find(x): 

if x^ x.parent then 

x.parent FIND {x.parent) 
return x.parent 



Figure 3: The union-find methods makeset, unite, and find with rank heuristic 
and path compression. Here, UNITE (x,};) returns X,F so that X is the new represen- 
tative of [x] U \y\ , and Y is the old representative of either [x] or \y\ which is no longer 
used. 

respective next and previous element in the doubly-linked outlist and inlist where e 
appears. In addition, e contains the links to two nodes e.tail and e.head, which never 
change, so that e is always an edge from e.tail to e.head in the original graph. In the 
current reduced graph at any stage of the computation, e is an edge from FIND (e. to//) 
to FlND{e.head), so these fields of e are not updated when e is moved to another 
node in an edgelist; this allows to move all incident edges from one node to another 
in constant time. 

Figure 4 gives pseudo-code for the contraction (a) described above. The three 
nodes U^V^W are the representatives of their partition classes, and only for these 
nodes the lists of outgoing and incoming edges and their matched edge are relevant. 
The unmatched edge e appears in U. outlist and has head V, so that U — find (^. to//) 
and V = FlNT>{e.head) = FlND(m.to//) , even though it may be possible that e.head ^ 
m.tail like for x,y in Figure 5. The matched edge m from V to W is obtained as 
V.matched (and equals W. matched), because V. outlist is empty so V has no outgoing 
unmatched edge (but has to have an outgoing edge due to the Eulerian orientation). 

After the shrink operation, the reduced graph no longer contains the edges e 
and m and the node V. (However, these are preserved for later re-insertion, helped 
by the field m.partner assigned to e in line 6 of SHRINK, discussed below along with 
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SHRlNK(e,m): 

1 U,W ■(r- FIND {e. tail), FIND {m. head) 

2 remove e from U.outlist 
*3 V -i^ FIND (e. head) 

* 4 sleepcounter sleepcounter + 1 

* 5 V.sleeptime sleepcounter 
6 m.partner ^ e 

1 Z,y^ UNITE (f/,VK) 

8 append Y.outlist to X.outlist 

9 append Y.inlist to X.inlist 

10 X. matched •<— U. matched 



Figure 4: The shrink operation that removes the unmatched edge e from U to V 
and matched edge m from V to W from the current graph and merges the edgeUsts 
of U and W. The code in the starred Unes 3-5 is only needed to reason about the 
method and can be omitted. 




Figure 5: The equivalence classes [U], [V], [W] and edges e and m in the shrink 
operation. A wiggly line denotes a matched edge. 

lines 3-5.) The edge e is removed from the list of outgoing edges of U in line 2. The 
equivalence classes for U and W are united in line 7 where either U or W becomes 
the new representative, stored in X . The lists of outgoing and incoming edges of the 
representative Y that is no longer in use are appended to those of X in lines 8 and 9. 
A node can only lose but never gain the status of being a representative, so there is 
no need to delete the edgelists of 7. If the new representative X is W, its current 
matched edge m has to be replaced by the matched edge U. matched as in line 10 
(which has no effect ifX = U). The Euler property of the reduced graph is preserved 
because the outdegree of X is the sum of the outdegrees of U and W minus one, and 
so is the indegree (the missing edges are e and m). 

The list operations in lines 2, 8, 9 of SHRINK can be performed in constant time. 
For that purpose, it is useful to store the lists of outgoing and incoming unmatched 
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m 



Edge fields: 



partner 




nextout 






prevout 


Z 




tail 


y 




head 






nextin 






previn 



Figure 6: Example of a graph with the out- and inlists for the nodes x^y^z accessed by 
sentinels (dummy edges) Sx,Sy,s^ shown in gray. They use the same fields nextout, 
prevout, nextin, previn as the unmatched edges ei,e2)^3 except for tail and head 
which are ignored. The matched edge m is stored directly, linked to by y.matched 
and z.matched, and not in a list. The m.nextout field can be used to link to m.partner. 



edges of a node u as doubly-linked circular lists that start with a "sentinel" (dummy 
edge), denoted by Su (see Gormen et al., 2001, Section 10.2). Figure 6 gives an 
example of small graph (which is neither Eulerian nor has a perfect matching). The 
three unmatched edges are ^1,^2,^3 and the matched edge is m. The outlist of x 
contains 61,^2. the outlist of y is empty, and the outlist of z contains ^3. The inlist 
of each node contains exactly one edge. The first and last element of the outlist of 
a node u is pointed to by Su-nextout and Su-prevout, which are both itself when 
the list is empty, as for m = 3; in the example. The inlist is similarly accessed via 
Su-nextin and Su-previn. Each append operation in line 8 or 9 of SHRINK is then 
performed by changing four pointers. The remove operation in line 2 can, in fact, be 
done directly from e, again by changing four pointers, here of the next and previous 
edge in the list (which may be a sentinel). Due to the sentinels, e does not need the 
information of which node it is currently attached to, so line 2 should be written (a 
bit more obscurely) as "remove e from its outlist" (that is, the outlist it is currently 
contained in), without reference to U . 

Figure 7 shows the whole algorithm that finds a perfect matching of opposite 
sign via a sign- switching cycle. Initialization takes place in lines 1-6, which will be 
explained when the respective fields and variables are used. 

The main computation starts at step A. The first node V is the head of a matched 
edge. This assures that, due to the Euler property, this node has at least one outgoing 
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FIND 


.OPPOSITELY_SIGNED_MATCHING : 


1 


for all nodes u 


2 


makeset(m) 


3 


u.origmatched •<— u.matched 


4 


u.visited 


* 5 


u.sleeptime 


* 6 


sleepcounter 


A : 


m any matched edge of current graph 


7 


V FmD{m.head) 


8 




B : 


visitednode[vc] •<— V 


9 


V. visited vc 


10 


if V.outlist is not empty then 


11 


e ^ first edge in V.outlist 


12 


W FIND(e./iearf) 


13 


visitededge[vc] •<— e 


14 


vc vc + 1 


15 


checkvisited(W) 


16 




17 


go to B 


18 


else 


19 


m ■<— V.matched 


20 


W ^ FlND{m.head) 


21 


vc vc — 1 


22 


?7, e •<— v/5/?e<inorfc[vc] , visitededge[vc] 


23 


it W = U then 


24 


return expandc ycle (e, m) 


25 


SHRINK(e,m) 


26 


checkvisited (ly) 


27 


11 vc > 1 then 


28 


V ^ nND(W) 


29 


go to B 


30 


else 


31 


go to A 



Figure 7: The main method FIND_0PP0SITELY_SIGNED_MATCHING for an Euler 
graph with a given perfect matching. 

unmatched edge that may be the first edge e of an edge pair c,m that is contracted 
with the SHRINK method. Starting from step B, a path of unmatched edges is grown 
with its nodes stored in visitednode[l] , . . . , visitednode[vc] where vc counts the num- 
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ber of visited nodes, and edges stored in visitededge[l] visitededge[vc — 1] , where 
visitededge[i] is the edge from visitednode[i] to visitednode[i + 1] for I <i <vc. A 
node u is recognized as visited on that path when u.visited is positive, which is the 
index i so that u = visitednode[i]. This field is initialized in line 4 as initially zero 
(unvisited). 

Line 10 tests if V has a non-empty list of outgoing unmatched edges, which is 
true when vc = 1 . The next node, following the first edge e of that list, is W . Line 15 
checks with the method CHECKVISITED, shown in Figure 8, if W has been visited 
before. If that is the case, then all edges in the corresponding cycle are completely 
removed from the graph and the nodes are marked as unvisited (lines 3 and 4 of 
CHECKVISITED in Figure 8), and vc is reset to the beginning of that cycle. In any 
case, W is the next node of the path, and the loop repeats at step B via line 17. 



CHECKVISITED (W): 

1 if W. visited > then 

2 for / ^ W.visited, . . . , vc — 1 

3 remove visitededge[i] from its outlist and inlist 

4 visitednode [i] . visited •<— 

5 vc •<— W.visited 



Figure 8: The CHECKVISITED method that checks if node W has already been vis- 
ited, and if yes deletes the encountered cycle of unmatched edges and updates vc. 

Lines 1 8-3 1 deal with the case that V has no outgoing unmatched edge, which can 
only hold if vc > 1 . Then the matched edge m incident to V is necessarily outgoing 
due to the Euler property and because V has an incoming edge e from UtoV, which 
is found in line 22. This edge is normally removed in the SHRINK operation and then 
no longer part of the path, which is why vc is decremented in line 21 (node V will no 
longer be part of the graph and can keep its visited field). However, a sign-switching 
cycle is found ifW — U (see Figure 9), which is tested in line 23 and dealt with in the 
EXPANDCYCLE method called in line 24, which terminates the algorithm and will be 
explained below. 

If W ^ U, then SHRlNK(c,m) is called in line 25. Afterwards, node W is still 
the old representative of the head node of m, as it was used in finding the path of 
unmatched edges. Node W may be part of that path, as tested (and the possible cycle 
removed) in line 26. If W has been visited, then W.visited < vc because W ^ U , 
so at least one edge is removed. If vc > 1 (which holds in particular if W has not 
been visited), then the path is now grown from FIND(V7) in line 28, where the FIND 
operation is needed to update visitednode[vc] in step B because the old representative 
U may have been changed to W after the unite operation in line 7 of SHRINK in 
Figure 4. 

The case W.visited = 1 needs special treatment, which results in vc = 1 and hap- 
pens when m is same as the initial matched edge m in step A. In that case, m is re- 
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Figure 9: The equivalence classes [?7], [V] when a sign-switching cycle has been 
found. 

moved via SHRINK, and FIND(W) may now be the tail rather than head of a matched 
edge, and then may no longer have an unmatched outgoing edge, which is necessary 
for lines 21-22 to work. For that reason, the loop goes back to A rather than B, as 
in line 31; note that W. visited has been set to zero in line 4 of CHECKVISITED with 
i = vc = I. In step A, a new matched edge that has not yet been removed in a call 
to SHRINK can be found in constant time by storing all matched edges in a doubly- 
linked list (for example, using the fields nextin and previn that so far are unused for 
matched edges, see Figure 6); a matched edge m should be deleted from that list after 
line 6 of shrink in Figure 4, for example. 



EXPANDC YCLE (e, m) : 



1 C-(-{(e,m)} 

2 RECONNECT (e./zea J, m. to//, C) 

3 RECONNECT (e. to//, m./ieaii,C) 

4 for all (e,m) e C 

5 make e a matched edge and m an unmatched edge 

6 return graph with this new matching 

RECONNECT (x,);, C) : 

7 if jc 7^ y then 

8 m •<— x.origmatched 

9 e m.partner 

10 C^CU{(£',m)} 

11 RECONNECT {e. head, m.tail,C) 

12 RECONNECT (e. to//, y,C) 



Figure 10: The expandcycle and the recursive RECONNECT method that create 
the sign-switching cycle and with it the oppositely signed matching. 
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We now discuss how to re-insert the contracted edges into the graph once a 

sign- switching cycle has been found, which is done in the expandcycle method 
in Figure 10. The method itself is straightforward. Recall that edges of the cur- 
rent graph are stored with the representatives of equivalence classes, where an un- 
matched edge is accessed in line 1 1 and a matched edge in line 19 of the main method 
FlND_OPPOSlTELY_SlGNED_MATCHlNG in Figure 7. The sign-switching cycle will 
be reconstructed using the original endpoints of the edges. For each node u, the orig- 
inal matched edge incident to u is stored in u.origmatched (see line 3 of the main 
method), because u.matched may be modified (in line 10 of the shrink method). 

In order to explain the expandcycle method, we record the time at which the 
SHRINK operation has been applied to a node V. This is done in lines 3-5 in Figure 4 
using using the field V.sleeptime and the global variable sleepcounter, which are 
initialized in lines 5-6 of Figure 7. These lines have a "star" to indicate that they do 
not affect the algorithm, and can therefore be omitted. We use them to reason about 
the correctness of the expandcycle method. 

The contraction SHRINK (e,m) affects three equivalence classes [U], [V], [W] with 
representatives ?7,y, W as shown in Figure 5. All nodes in [V] become inaccessible 
afterwards, but the equivalence class still exists (and is in fact still represented in 
the union tree by those nodes v so that V = find(v), although the union-find data 
structure will no longer be used for these nodes). We say that all nodes in [V] become 
asleep at the time recorded in the positive integer V.sleeptime. Any node u so that 
FlNT){u) .sleeptime = is called awake. 

Lemma 15 During the main method find_oppositely_signed_matching, the 
following condition holds after any statement from step A onwards. Let [U] be an 
equivalence class of nodes with representative U and let m' = U.matched. Then 
there is exactly one node u in \JJ] and another node z not in [U] so that: 

(i) IfU is awake, then z is awake and {m,z} = {m' .head, m' .tail}. 

(ii) IfU is asleep, then z is awake or asleep with later sleeptime than U , and u — 
m' .tail, z — m' .head. 

In either case: 

(Hi) For every node y in \U] — {a}, let m — y.origmatched. Then y = m.head, the 
node m.tail is asleep (with earlier sleeptime than U if U is asleep), there is an 
edge e so that e — m.partner, the nodes m.tail and e.head are equivalent and 
with X = e.tail we have xE [U] — {y}. 

Proof. We prove this by induction over the number of calls to the shrink method, 
which are the only times when the equivalence classes change. Initially, all equiv- 
alence classes are singletons and all nodes are awake. Then \U\ = {U}, only (i) 
applies, where m' is the matched edge incident to u which is either the tail or head of 
m' and z is the other endpoint of m' , and (iii) holds trivially. 

Figure 1 1 shows the general case of the lemma. Consider the three equivalence 
classes [U ] , [V] , [W] with representatives U, V, W as shown in Figure 5 in the notation 
used for the shrink method. Before shrink (e,m) is called, the lemma applies by 
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Figure 11: Illustration of Lemma 15. The matched edge m' with endpoints u and z 
may have either orientation if u is awake, otherwise u = m' .tail as stated in (ii). 

inductive assumption to each of the three classes \U], [V], [W] in place of [[/]. The 
unique matched edge m' that goes outside the equivalence class to an awake node 
is m for [V ] and [V7] , and for \U ] it is some other matched edge m' (not shown in 
Figure 5) which will be that edge after [U ] and [VF] have been united. There is no 
edge other than e or m from a node in [V] to an awake node outside [V] because only 
in this case (when V has in- and outdegree one in the reduced graph) the SHRINK 
method is called. Every node in \U], [V] , or [VF] is the endpoint of a matched edge in 
the original graph, and other than the endpoints of m and m' they are all equal to the 
head of such a matched edge, with its tail node asleep, by inductive assumption (iii). 

After the shrink operation, [U] and [W] become a single equivalence class, and 
all nodes in [V] becomes asleep. The only node y in the new class \U] U [V7] for 
which (iii) does not hold by inductive assumption is m.head, but then e takes exactly 
the described role as m.partner. In particular, x — e.tail ^ y, because x ^[U] and 
y e [VF] and [U ] ^ ^] . In addition, [V ] changes its status from awake to asleep, and 
all nodes in [V] — {m.tail} are heads of matched edges that connect to equivalence 
classes that went asleep before V as claimed in (iii) by the inductive hypothesis. This 
completes the induction. □ 

The previous lemma implies that any two endpoints of a matched edge belong 
to different equivalence classes. A key observation in (iii) is that for any y in an 
equivalence class [U] that is not the endpoint u of the "awake" matched edge m' 
there is another node x different from y in that class (which may be u) given by 
X = y.origmatched.partner.tail. 

Lemma 16 Consider nodes x,y and a set C with the following properties: x and y 
are equivalent, and x is the endpoint of an unmatched edge and y is the endpoint 
of an oppositely oriented matched edge taken from the pairs of unmatched-matched 
edge pairs in C, as in (i) or (ii) in Figure 12. Then after RECONNECT (x,3;,C), the 
new edges in C form a path of alternating matched-unmatched edges from xtoy with 
the same number of matched and unmatched forward-pointing edges. 
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Proof. If x = y, then the claim holds trivially with a path of zero length because no 
edge pairs are added to C. Otherwise, we apply Lemma 15 to the equivalence class 
that contains x and y, as shown in the right picture in Figure 12, where x takes the 
role of y in Lemma 15(iii). Hence, node x is the head of some matched edge m with 
partner e, and x' = e.tail ^ x. Then in the method RECONNECT (line 10 in Figure 10), 
(e, m) is added to C. We then apply the claim recursively to e.head, m.tail instead of 
and then Xo x' ,y instead of x,y, where the assumptions apply, exactly as in lines 
1 1 and 12 of reconnect. So there are alternating paths as described from e.head to 
m.tail and from x' io y. The resulting path from x to y composed of these paths and 
the edges m and e has the same number of forward-pointing matched and unmatched 
edges, because m and e point in the same direction (in this case, backwards) along 
the path. □ 

In the EXPANDCYCLE method, lines 2 and 3 in Figure 10 call reconnect (x, y, C) 
for the endpoints x, y of the unmatched-edge pair shown in Figure 9 that results when 
EXPANDCYCLE is called from the main method (line 24 of Figure 7), first for x,y in 
[y] and then for x,y in [U] in Figure 9. In both cases, Lemma 16 applies, and the 
paths together with the first edge pair (e,m) form a sign-switching cycle. 

Finally, exchanging the matched and unmatched edges as in line 5 of EXPAND- 
CYCLE can be done as described, irrespective of the order of the edges in the cycle, 
just using the pairs (e,m) in C (which are oriented in the same direction along the 
cycle), which suffices to obtain a matching of opposite sign. 

This concludes the detailed description of the algorithm. It has near-linear run- 
ning time in the number of edges, because each unmatched edge is visited at most 
once and either discarded or contracted in the course of the algorithm. 

We illustrate the computation with an example shown in Figure 13. Suppose that 
edge lists contain edges in alphabetical order. The first node is 1. The first three 
iterations follow the unmatched edges a,b,c, so that vc and the arrays visitednode 
and visitededge have the following contents: 
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Figure 13: Example to illustrate the algorithm. Unmatched edges are marked ato h, 
matched edges are identified by their endpoints. The first matched edge is 61 . 

vc = 4, visitednode = [\ 2 3 4], visitededge = [a b c]. 

Node 4 has an empty outlist, so that the computation continues at line 18 (all line 
numbers refer to the main method FlND_OPPOSlTELY_SlGNED_MATCHlNG in Fig- 
ure 7 unless specified otherwise). The matched edge is m = 45 with endpoint W = 5, 
and SHRINK (c, 45) is called in line 25. Afterwards, 

A5 .partner = c , vc = 3, visitednode = [\ 1 3], visitededge = [a b] 

and the reduced graph is shown in Figure 14. The nodes 3 and 5 have been united into 
the equivalence class {3,5} which has representative 5 because the unite operation 
in Figure 3 chooses the second representative Y if the original representatives have 
equal rank. The outgoing edges from node 5 are J and h in that order, because the 
outlist of 3 has been appended to that of 5. 

In line 26, checkvisited(W) has no effect because W. visited = 0. In line 28, 
V FiND(iy) = = 5, so that after going back to step B, 

vc = 3, visitednode = [I 2 5], visitededge = [a b]. 

Line 1 1 then follows the unmatched edge d with 

vc = 4, visitednode = [I 2 5 6], visitededge = [a b d] 

after which again a contraction is needed, because V.outlist is empty, this time with 
U,V,W = 5,6,1 and call to SHRINk(J,61). In the shrink method, UNITe(5,1) 
returns 5, 1 because 5. rank = 1 > = I. rank. The resulting graph, after line 25 is 
completed, is shown on the left in Figure 15, where 

6l.partner = d, W = l, vc = 3, visitednode = [I 2 5], visitededge = [a b]. 
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Figure 14: The reduced graph after the first contraction with shrink (c, 45). The 
equivalence class with nodes 5, 3 is written with the representative listed first. 

Now consider the call to CHECKVISITED (W) in line 26 and note that W is still the old 
node 1 used before the contraction; recall that this is done because that representative 
is possibly stored in the visitednode array, which it indeed is at index W. visited — 1 . 
The deletion of the detected cycle of unmatched edges in lines 3 and 4 of CHECKVIS- 
ITED (see Figure 8) then produces the reduced graph shown on the right in Figure 15. 




Figure 15: Left: after SHRiNK(ci,61), right: after checkvisited(I). 



Normally, the next node V would be find(W) in line 28. However, the case 
vc = 1 applies (recall the reason that the incoming matched edge of visitednode[l] 
has been removed by the contraction), and so the computation continues via line 3 1 
to step A. Suppose the first node is now 2. Then the computation follows edges e,f, h 
and reaches node 8, with 

vc = 4, visitednode =[2 7 5 8], visitededge =[e f h]. 

Contraction with SHRINK(/j, 87) gives the graph shown on the left in Figure 16 where 

S7 .partner = h, vc = 3, visitednode = [2 1 5], visitededge = [e f] 
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and after checkvisited (7) the edge / is removed, resulting in 

vc = 2, visitednode = [2 7] , visitededge = [e] . 

This time, vc > 1, so with V ■<— find(W) = 5 we go back via Une 29 to step B, after 
which 

vc = 2, visitednode = [2 5] , visitededge = [e] 
with the graph on the right in Figure 16. 




Figure 16: Left: reduced graph after the third contraction shrink(/i, 87) , right: after 

CHECKVISITED (7). 

Now V has only the outgoing unmatched edge g, giving in line 14 (where the last 
entry visitednode[vc] has not yet been assigned) 

vc = 3, visitednode = [2 5], visitededge = [e g] 

and removal of the edge g gives the graph on the left of Figure 17. At the next 
iteration, after line 9, 

vc = 2, visitednode = [2 5] , visitededge = [e] 

where S.outlist is empty, m = 52, and now W = 2 = U = visitednode[\] in line 23. 
Now the final stage of the algorithm is called in line 24 with EXPANDCYCLE (e, m) . 
The original endpoints of the two edges are (see Figure 13): e.tail = 2, e.head = 7, 
m.tail — 3, m.head — 2. Line 2 of Figure 10 makes the call RECONNECT (7, 3, {e,m}) 
which is nontrivial because with j:, j = 7 , 3 we have x^yin line 7 of Figure 10. With 
x.origmatched = 87 and S7. partner = h, we get C = {(e, 52), (h, 87)} (where 52 is 
just our current name for the matched edge, with its original endpoints it is the edge 
32). All other calls to RECONNECT (;c,);,C) then have no effect because x = y. The 
resulting sign- switching cycle C is shown on the right in Figure 17. 

In this example, every edge of the graph is visited during the algorithm, and the 
reduced graph at the end consists justs of the oppositely oriented unmatched and 
matched edge that define a trivial sign-switching cycle. The original graph in Fig- 
ure 14 already has such an edge pair in the form of b, 32, which is not discovered by 
the described run of the algorithm. Here, not all matched edges and their partners that 
have been removed by the shrink operation are used (namely not 45, c and 61, J). 
In other cases, the algorithm may also terminate with parts of the graph left unvisited, 
or unmatched edges in the visitededge array that are not removed. 
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Figure 17: Left: graph after the removal of g, which has a sign-switching cycle. 
Right: cycle after the calls to reconnect in expandcycle. 
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